UNCLASSIFIED 


ft  PROBABILISTIC  HODEl  OF  THE  APPARENT  RADIANCE  OF  A 
ROUGH  SEA(U)  NAVAL  RESEARCH  LAB  HASHINGTON  DC 
R  0  PRIEST  ET  AL.  S4  AAR  88  NRL-HR-8892 


aval  Research  Laboratory 


Washington,  DC  20375-5000 


NRL  Memorandum  Report 


Js ilL  •  i'.L 


A  Probabilistic  Model  of  the  Apparent 
Radiance  of  a  Rough  Sea 

R.  G.  Priest  and  I.  B.  Schwartz 

Advanced  Concepts  Branch 
Optical  Sciences  Division 


rypr* 

£LECTE 

APR  2  5  1988 


March  4,  1988 


Approved  for  public  release;  distribution  unlimited 


V 


REPORT  DOCUMENTATION  PAGE 


lb  RESTRICTIVE  MARKINGS 


form  Approved 
OMB  No.  0704-0! SB 


la.  REPORT  SECURITY  CLASSIFICATION 

UNCLASSIFIED 

2a.  SECURITY  CLASSIFICATION  AUTHORITY 

Zb.  DECLASSIFICATION  /  DOWNGRADING  SCHEOULE 

4.  PERFORMING  ORGANIZATION  REPORT  NUMSEr7s>  ~| 

NRL  Memorandum  Report  6092 

6a.  NAME  OF  PERFORMING  ORGANIZATION 

Naval  Research  Laboratory 

6b  OFFICE  SYMBOL 
(If  applicable) 

Code  6522 

6c  AOORESS  (Oty,  State,  end  ZIP  Code) 

Washington,  DC  20375-5000 

8a.  NAME  OF  FUNOING/ SPONSORING 
ORGANIZATION 

Office  of  Naval  Technology 

8b  OFFICE  SYMBOL 
(If  applicable) 

ONT 

8c  AOORESS  (City,  State,  end  ZIP  Code) 

Arlington,  VA  22217-5000 

3  DISTRIBUTION /AVAILABILITY  OF  REPORT 

Approved  for  public  release;  distribution 
unlimited. 


5  MONITORING  ORGANIZATION  REPORT  NUMBER(S) 


7*.  NAME  OF  MONITORING  ORGANIZATION 

Office  of  Naval  Technology 


7b  AODRESS  (Cty,  Sun.  too  ZIP  Cod*) 

Arlington,  VA  22217-5000 


9  PROCUREMENT  INSTRUMENT  IDENTIFICATION  NUM8ER 


10  SOURCE  OF  FUNDING  NUMBERS 


PROGRAM 
ELEMENT  NO 

62712N 


PROiECT 

TASK 

NO 

NO 

WORK  UNIT 
ACCESSION  NO 

65139 


1 1 .  TITLE  (Mud*  Socunty  Claaification) 

A  Probabilistic  Model  of  the  Apparent  Radiance  of  a  Rough  Sea 


12.  PERSONAL  AUTHOR(S) 

Priest,  R.G.  and  Schwartz,  l.B. 


13a.  type  OF  REPORT  113b  TIME  COVERED 

Memo  I  from  10/86  to  10/87 


14  OATE  OF  REPORT  (Year,  Month.  Day)  15  PAGE  COUNT 

1988  March  U  16 


COSATI  CODES 


SUB-GROUP 


18  SUBJECT  TERMS  (Cootmut  on  reverse  if  necesSJry  end  idtntify  by  Mock  number) 

Infrared  Rough  surface 

Probability  Rough  sea 


19  ABSTRACT  (Continue  on  reverse  if  necessary  and  identify  by  block  number) 


A  widely  used  approach  to  the  specification  of  the  background  statistics  of  wind  driven  sea  states  is  the 
power  spectral  method  (*).  Coupled  with  a  ray  trace  of  the  sky  radiance,  one  can  create  deterministic  images 
of  a  rough  sea  in  the  infrared.  In  this  paper  a  probabilistic  model  of  a  rough  sea  is  presented  as  an  alterna¬ 
tive  to  the  power  spectral  approach.  The  model  takes  into  account  self-shadowing  of  the  sea  surface  which  is 
important  when  the  sea  is  viewed  at  near  grazing  angles.  Given  such  a  model,  the  apparent  radiance  of  a 
rough  sea  is  shown  to  emit  a  substantially  lower  amount  of  radiation  when  compared  to  a  smooth  sea.  The 
sea,  being  a  poor  emitter  of  radiation  at  near  grazing  angles,  reflects  cool  sky  radiation  from  higher  points  in 
the  sky  and  therefore  looks  colder  than  its  blackbody  temperature. 

(•)!.  B.  Schwartz  and  D.  Han,  Emiuivity  as  a  function  of  surface  roughness.  A  computer  model.  NRL  Memorandum  Report  No. 
3816.  1986. 


20  DISTRIBUTION  AVAILABILITY  OF  ABSTRACT  21  ABSTRACT  SECURITY  CLASSIFICATION 

bduNCLASSlF'ED/UNLiMITED  □  SAME  AS  PPr  □  3' C  j$E»S  UNCLASSIFIED  _ 


22a  NAME  OF  RESPONSIBLE  'NOIVIOUAL 

I.  B.  Schwartz 


DO  Form  1473,  JUN  86  Previous  tuitions  are  obsolete  _ ^ec,°itv  c.asSiFCatqn  Of  -aq; 

S/M  0102-LF-0U-6603 


mnxm-p 


ScC^Pl jy  ClASS. P  CAT  on  Qf 


ws 


raw 


rtr. 


TO 


\ 


9 


m 


alk  *  - 


1'*  I  ,  t 


i  la*  U  U  >1  *»<  J»*J 


rA  ■**  *ith  «*fc 


A  PROBABILISTIC  MODEL  OF  THE 
APPARENT  RADIANCE  OF  A  ROUGH  SEA 


INTRODUCTION 


The  infrared  signature  of  ships  is  most  often  described  in  terms  of  contrast  to  the  background 
infrared  signature  of  the  sea.  Since  the  contrast  signature  involves  the  signature  of  the  sea,  an  impor¬ 
tant  task  for  signature  modeling  is  to  calculate  the  radiance  of  the  sea  surface  as  a  function  of  viewing 
angle.  This  is  an  interesting  problem  because  of  the  complication  caused  by  the  roughness  of  the  sea 
surface.  Since  the  infrared  wavelength  is  much  shorter  than  the  length  scale  of  disturbances  on  the 
water’s  surface,  the  rough  sea  may  be  thought  of  as  an  ensemble  of  small  facets  each  of  which  emits 
and  reflects  (from  the  sky,  and  possibly  other  facets)  infrared  radiation.  In  cases  where  the  length 
corresponding  to  the  projection  of  the  instantaneous  field  of  view  (IFOV)  of  the  sensor  onto  the  ocean 
surface  is  greater  than  the  length  scale  of  sea  disturbances,  reflection  and  emission  properties  may  be 
calculated  by  taking  an  average  over  facets  weighted,  in  an  appropriate  manner,  by  the  probability 
density  function  for  facet  orientation. 

It  is  easy  to  see  that  the  distribution  of  facet  orientations  can  have  a  large  effect  on  the  infrared 
radiation  emanating  from  the  sea  surface.  If  the  sea  were  smooth,  all  of  the  facets  would  be  horizon¬ 
tally  oriented.  However  when  a  rough  sea  surface  is  viewed  near  grazing,  most  of  the  facets  seen  by 
the  sensor  are  tipped  away  from  the  horizontal  toward  the  sensor  [1,2].  This  phenomenon  impacts 
upon  both  the  amount  of  radiation  emitted  from  the  sea  and  the  amount  of  radiation  reflected  from  the 
sky.  The  average  dlt  angle  causes  radiation  from  elevation  angles  relatively  higher  in  the  sky  to  be 
reflected  than  would  be  the  case  if  the  sea  where  not  rough.  Since  the  strength  of  sky  radiation 
decreases  sharply  with  increasing  elevation  angle  the  effect  of  sea  roughness  is  to  decrease  the  amount 
of  radiation  reflected  from  the  sky.  At  the  same  time  the  decreased  angle  of  incidence  to  the  facet 
increases  the  emissivity  of  the  facet  [3].  This  effect  increases  the  amount  of  radiation  c mined  from 
the  facet.  At  angles  very  near  grazing  the  reflected  radiation  decrease  is  larger  than  the  emitted  radia¬ 
tion  increase.  For  other  angles  the  situation  is  reversed. 

A  major  theoretical  problem  is  that  some  facets  are  hidden  from  view  by  others.  This  is  called 
the  hidden  surface  problem  or  shadowing  problem.  Two  approaches  to  this  statistical  problem  are 
possible.  The  first  is  to  make  a  calculation  of  the  bidirectional  reflectance  distribution  function 
(BRDF)  under  the  assumption  that  there  is  no  shadowing  and  then  to  introduce  some  heuristic  correc¬ 
tions  to  take  shadowing  into  account.  This  approach  has  been  studied  by  Torrance  and  Sparrow  [4] 
for  the  case  of  rough  material  surfaces  and  by  Maxwell  [5]  for  the  case  of  rough  water.  The  second 
approach  involves  a  calculation  of  the  contribution  to  the  total  surface  area  seen  by  a  sensor  that  is 
made  by  facets  of  a  given  orientation.  This  approach  permits  a  simple  adjustment  to  be  made  for  one 
of  the  effects  of  shadowing.  It  is  this  second  approach  which  is  presented  in  this  paper.  It  will  be 
shown  that  both  approaches  lead  to  the  same  result  in  cases  where  there  is  no  shadowing.  It  will  be 
seen  that  the  simple  correction  allowed  by  the  second  approach  leads  to  reasonable  results  for  viewing 
angles  near  grazing  incidence.  For  this  reason  it  represents  an  advance  over  the  corrected  BRDF 
approach. 
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STATISTICS  OF  FACET  PROJECTED  AREAS 

As  the  first  step  it  is  necessary  to  calculate  the  wave  slope  probability  density  function.  An  ana¬ 
lytic  result  is  available  for  the  case  in  which  the  wave  amplitudes  are  sufficiently  small  such  that 
linear  superposition  of  amplitudes  is  valid.  The  essential  elements  of  this  calculation  have  appeared  in 
the  literature  [6].  For  completeness  the  derivation  is  presented  in  appendix  A.  The  central  result  can 
be  expressed  in  terms  of  a  probability  density  for  the  orientation  of  a  unit  vector  n  normal  to  the  sur¬ 
face  of  a  facet.  The  simplest  description  of  this  density  is  in  terms  of  a  system  of  spherical  polar 
coordinates  with  polar  axis  (the  z  axis)  in  the  zenith  direction  and  azimuths  measured  with  respect  to 
the  upwind  direction  (the  x  axis).  If  0  and  are,  respectively  the  polar  and  azimuthal  angles  of  h , 
the  probability  density  is: 

-  (l/2T)(ab)-‘  tan  (0)  sec2  (0)*-“^ (l) 


-  j 


cos2  (<p)  sin2  (<fi) 


where  a2  is  the  variance  of  the  derivative  of  the  waveslope  in  the  upwind  direction  and  b-  is  the  vari¬ 
ance  in  the  crosswind  direction.  This  expression  may  be  thought  of  as  the  angular  probability  density 
for  finding  the  facet  at  a  given  place  on  the  surface  to  be  oriented  such  that  its  normal  has  the  angular 
coordinates  d,<p.  We  assume  there  is  no  correlation  between  upwind  and  crosswind  directions. 

The  basis  of  the  approach  developed  in  this  paper  is  the  calculation  of  a  fraction,  /(0,*>),  of  the 
total  projected  area  seen  by  a  sensor.  This  is  the  fraction  contributed  by  facets  falling  in  an  infini¬ 
tesimal  angular  range  centered  at  9,<p.  Once  this  fraction  is  determined  it  can  be  used  to  weight  the 
reflection  and  emission  contributions  of  the  facets  to  the  net  observed  radiance.  To  make  progress  it 
is  necessary  to  define  two  types  of  projected  areas.  The  first  type  is  horizontally  projected  area. 
This  is  the  horizontal  (parallel  to  the  mean  sea  level  )  area  obtained  by  projecting  a  facet  or  group  of 
facets  onto  the  horizontal  plane.  The  second  is  line-of-sight  (LOS)  projected  area.  LOS  projected 
area  is  the  area  obtained  by  projecting  a  facet  or  group  of  facets  onto  a  plane  normal  to  the  LOS  to 
the  sensor.  The  probability  density  function  p(9,<p)  is  related  to  the  horizontally  projected  area  in  the 
following  way:  the  horizontally  projected  area,  AA*,  contributed  by  facets  with  orientations  in  the 
infinitesimal  angular  range  A0A<*>  centered  at  0,<p  is  given  by  : 

AA*  -  A*p(0,*)A0A*  (2) 

where  A*  is  the  total  horizontally  projected  area  in  the  region  under  consideration.  Note  that  the 
integral  over  0  and  v>  of  the  RHS  of  Eq.  (2)  correctly  yields  Ah  due  to  the  fact  that  p{6.<p)  is  a  (prop¬ 
erly  normalized)  probability  density.  The  actual  (unprojected)  surface  area,  AA  ,  of  the  same  group  of 
facets  is  AA*  /cos  (0).  That  A A  is  greater  than  AA*  is  a  reflection  of  the  well  known  fact  that  the 
total  surface  area  of  a  rough  surface  is  greater  than  the  nominal  area  which  it  covers.  The  LOS  pro¬ 
jected  area  of  this  same  group  of  facets,  AALos>  ‘s  given  by  AA  cos  (fl),  were  Q  is  the  angle  between 
the  normal  to  the  facet  and  the  LOS  from  the  facet  to  the  sensor.  The  expression  for  A A  LOS  in  terms 
of  the  probability  density  function  is  thus: 

AA los  =  Ai,p(0,ifi)  - ~-A0A <p.  (3) 

cos  (0) 

The  total  LOS  projected  area  can  be  obtained  by  integrating  the  RHS  of  Eq.  (3)  over  0  and 
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Caution  must  be  exercised  in  the  use  of  this  result.  While  0  is  constrained  by  p  to  be  less  than 
t/2,  there  is  no  similar  constraint  on  0.  Consequently  the  RHS  of  Equation  (3)  can  be  less  than  zero 
for  some  orientations  of  the  LOS  to  the  sensor.  Eq.  (3)  associates  negative  projected  areas  with 
facets  that  are  oriented  such  that  their  normal  vector  points  away  from  the  sensor.  When  the  RHS  of 
Eq.  (3)  is  integrated,  the  negative  projected  areas  contributed  by  facets  that  point  away  from  the  sen¬ 
sor  partially  cancel  the  positive  projected  areas  of  other  facets. 

The  preceding  discussion  of  projected  area  makes  it  possible  to  give  an  expression  for  the  frac¬ 
tional  projected  area  / (0,?).  The  quantity  f(6,<fi)  ABA<fi  is  defined  to  the  fraction  of  the  LOS  pro¬ 
jected  area  seen  by  the  sensor  that  is  contributed  by  facets  with  orientations  in  the  range  A0A <p  cen¬ 
tered  at  8,<fi.  In  the  absence  of  shadowing  effects  /  is  given  by: 


P(9,<?) 


cos  (Q) 


dddtp 


/<••*> - ^rir — •  (4) 

j p(9,«>)  dOd* 

r  cos  (9) 

That  is,  the  expression  of  Eq.  (3)  is  divided  by  the  whole  projected  area  seen  by  the  sensor.  The  lim¬ 
its  of  the  integration  here  and  throughout  the  paper  are  understood  to  be  0  to  2x  for  the  azimuthal 
angles  and  0  to  r  for  the  polar  angles.  This  expression  must  be  modified  to  account  properly  for 
shadowing  effects.  The  form  of  this  expression  and  the  cautionary  observation  discussed  above 
together  suggest  a  heuristic  medication  to  Eq.  (4)  that  can  account  for  some  of  the  shadowing  effects. 
The  class  of  hidden  facets  of  interest  are  the  ones  that  are  turned  away  from  the  sensor.  These  facets 
arc  aways  hidden  from  view.  According  to  the  formulas  of  Eq.  (3)  and  (4)  these  facets  make  nega¬ 
tive  contributions  toward  projected  areas.  The  correct  contribution  for  these  facets  is  zero.  This 
observation  suggests  the  following  modified  form  of  Eq.  (4)  : 


/(«.*>  - 


P0*)  -  0) 

cos  (0) 


[ P(9,<p)  Hilt/ 2  -  awed* 

r  cos  id) 


Here  f  is  the  corrected  form  of  /  and  H  is  the  step  function  which  is  equal  to  unity  if  its  argument  is 
positive  and  zero  otherwise. 

An  alternative  method  for  arriving  at  Eq.  (3)  involves  consideration  of  a  shadowing  factor.  This 
is  the  probability  that  a  facet  of  given  orientation  is  not  hidden  from  view  of  the  sensor.  This  prob¬ 
lem  has  been  analyzed  in  the  context  of  RADAR  detection  [7].  The  result,  ignoring  correlation  of 
waveheight  and  slope,  is  that  the  shadowing  factor  is  proportional  to  Hir/2  -  Q).  The  proportional¬ 
ity  factor  is  independent  of  facet  orientation.  The  step  from  Eq.  (4)  to  Eq.  (5)  can  be  viewed  as  the 
replacement  of  LOS  projected  areas  with  the  product  of  LOS  projected  area  and  the  appropriate 
shadowing  factor. 

The  central  idea  introduced  in  this  paper  is  to  use  the  expression  of  Eq.  (5)  to  weight  facet  con¬ 
tributions  to  the  net  radiance.  For  purposes  of  comparison  with  the  results  of  the  BRDF  approach  it 
is  useful  to  note  that  the  integral  in  the  denominator  of  Eq.  (5)  can  be  explicitly  carried  out  in  the 
nonshadowing  case.  In  this  case  the  function  H  is  identically  equal  to  unity.  The  spherical  addition 
theorem  may  be  used  to  calculated  cos  (fi): 

cos  (0)  *  cos  (9)  cos  (0)  +  sin  (0)  sin  (0)  cos  (<p  -  £)  (6) 

where  0  and  £  are  the  polar  and  azimuthal  angles  of  a  unit  vector,  r,  which  is  parallel  to  the  LOS  to 


the  sensor.  (See  Fig.  1  for  specification  of  the  coordinate  system.)  The  integral  in  the  denominator  of 
Eq.  (5)  thus  involves  two  terms.  The  first  term  evaluates  immediately  to  cos  (0)  because  the  proba¬ 
bility  density  is  normalized.  The  second  term  involves  the  integral  of  the  product  of  cos  (<?  -  <p)  and 
p(8,<p).  The  probability  functions  considered  in  this  paper  all  have  inversion  symmetry.  That  is 
p(8,<p )  ■  p(8, r  -  <p).  With  this  condition  it  is  easy  to  see  that  the  second  integral  is  zero.  How¬ 
ever  it  can  be  shown  that  this  integral  vanishes  under  more  general  conditions.  The  integral  in  the 
denominator  of  Eq.  (5)  thus  evaluates  to  cos  8  in  the  nonshadowing  case.  Thus  in  this  case: 


f(8,<p)  -  /(0,*)  -  p(8,<p) 


cos  (0) 


cos(0)  cos  (0) 


(7) 


EMISSIVITY  AND  REFLECTIVITY 

The  facet  emissivity,  e,  and  reflectivity,  p,  are  both  functions  of  the  angle  Q.  The  Fresnel  for¬ 
mula  for  unpolarized  radiation  may  be  used  to  calculate  p(Q)  if  the  complex  optical  constant  for  water 
is  specified.  The  appropriate  value  is  obtained  by  computing  an  average  weighted  by  the  response 
function  of  the  sensor.  Kirchhoffs  law  in  the  form  e  =*  1  -  p  may  be  used  to  obtain  the  angularly 
dependent  emissivity.  The  intrinsic  radiance,  A,,  is  given  by  Af  *  e(Q)Afl  where  Rg  is  the  radiance 
of  a  blackbody  at  the  sea  temperature.  The  reflected  radiance,  Rr  is  given  by  Rr  =  p(Q)A,(0'), 
where  R,  is  the  radiance  downwelling  from  the  sky  (sky  radiance).  The  sky  radiance  depends  on  the 
angle,  8',  between  the  unit  vector  r'  and  the  zenith  direction.  The  unit  vector  r'  is  parallel  to  the  ray 
generated  by  the  reflection  of  the  LOS  from  the  sensor  by  the  facet.  It  has  the  angular  coordinates  d' 
and  <p'.  A  vector  expression  can  be  given  for  r'  in  terms  of  n  and  r: 

?'  —  2  cos  (0 )n  -  r.  (8) 

An  expression  for  8'  in  terms  of  Q  and  6  can  be  obtained  from  this  equation  by  taking  the  inner 
product  of  both  sides  of  the  equation  with  the  unit  vector  pointing  in  the  zenith  direction: 

cos  (8')  =  2  cos  (Q)  cos  (0)  -  cos  (0).  (9) 

ROUGH  SEA  RADIANCE 

The  apparent  radiance  of  the  rough  sea,  Ra,  is  the  expected  value  of  the  sum  of  R,  and  Rr . 
This  expected  value  is  given  by: 

Aa(0)«  \d8\dpf  (A,  +  Rr) 

=  \d8\d<pf(8,p)\[l  -  p(Q))RB  +  p(Q)A,(0')i.  (10) 

Here  the  form  of  /  is  given  by  Eqs.  (5)  and  (1),  1)  is  given  by  Eq.  (6)  and  0'  is  given  by  Eq.  (9). 
The  Fresnel  formula  is  to  be  used  for  the  functional  form  of  p(Q)  and  the  functional  form  of  Rs  must 
be  externally  specified.  An  atmospheric^  modeling  code  such  as  Lowtran  VI  can  be  used  to  create  a 
tabulated  form  of  this  function.  Since  /  is  a  normalized  weighting  function,  it  can  be  verified  that 
Eq.  (10)  satisfies  an  important  consistency  relation.  Specifically,  in  the  case  where  R,  is  a  constant 
equal  to  RB,  one  finds  that  Ra  *  Ra.  This  result  is  required  by  the  second  law  of  thermodynamics. 
It  forms  a  useful  check  on  the  consistency  of  the  general  approach. 
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As  mentioned  in  the  introduction  it  can  be  shown  that  the  results  obtained  in  this  paper  agree 
with  the  results  obtained  using  the  BRDF  approach  in  the  case  where  there  is  no  shadowing.  In  this 
case  /  is  given  by  Eq.  (7).  If  this  form  is  substituted  into  Eq.  (10)  with  the  p{0,<p)  given  by  Eq.  (1) 
with  a2  *  b2  m  a2/!  (isotropic  case)  the  weighting  terms  in  Eq.  (10)  become: 

tan  (0)  cos  (0)  exp  (-tan2  (6)/ a2) 
cos3  (0)  cos  (0) 

Apart  from  notational  differences,  this  weighting  factor  is  the  same  as  the  one  in  [8]. 

RESULTS 

Equation  (10)  is  the  central  result  of  this  paper.  To  see  the  typical  magnitude  of  the  effects 
caused  by  a  rough  sea,  a  functional  form  for  R,  must  be  chosen.  As  an  example  we  have  exercised 
the  Lowtran  program  with  summer  midlatitude  parameters  and  have  obtained  the  long  wave  IR  (8-12 
micron)  sky  radiance  form  shown  in  Fig.  2.  For  convenience,  the  value  of  the  blackbody  radiance  of 
the  sea,  RB,  which  appears  in  Eq.  (10)  was  taken  to  be  equal  to  the  value  of  Rt  at  the  horizon.  This 
corresponds  to  selecting  temperatures  of  the  sea  and  of  the  air  just  above  the  sea  to  be  equal.  These 
Lowtran  results,  together  with  Eq.  (10)  (with  a2  =*  b2  m  a2/!)  were  used  to  calculate  Z(  0')  *  1  - 
Ra(0')/Rg  for  a  set  of  values  of  a.  For  this  calculation  the  optical  constants  appropriate  to  the  long 
wave  IR  were  used  in  the  Fresnel  reflectivity  formula.  Since  there  is  an  empirical  relationship 
between  a  and  the  wind  speed  given  by  [9]: 

a2  -  .003  +  .00512  w  (12) 

where  w  is  the  wind  speed  in  m/sec,  this  set  of  results  is  effectively  parameterized  by  the  wind  speed. 
Figure  3  presents  the  calculated  function  Z(0')  for  selected  values  of  cr.  It  can  be  observed  that  for 
values  of  O'  near  75°  (15°  away  from  grazing)  modest  values  of  the  wind  speed  produce  substantial 
effects.  The  sea  can  appear  to  emit  up  to  10%  less  radiation  than  a  calm  sea  would  appear  to  emit. 
Effects  of  this  magnitude  can  have  a  substantial  impact  on  ship  detectability. 

DISCUSSION 

In  this  paper  we  have  presented  a  probabilistic  model  of  the  apparent  radiance  of  a  rough  sea 
surface.  It  is  observed  that  when  compared  to  a  smooth  sea  surface,  a  rough  sea  can  emit  a  substan¬ 
tially  lower  amount  of  radiation  in  the  infrared.  The  sea,  being  a  poor  emitter  of  radiation  at  near 
grazing  angles,  reflects  cool  sky  radiation  from  higher  points  in  the  sky  and  therefore  looks  colder 
than  its  blackbody  temperature.  Other  authors  have  modelled  this  phenomenon  as  well  but  without 
the  aid  of  a  probability  model  of  the  sea  slopes  [1,2].  In  some  extreme  cases,  other  authors  rely  on 
one  emissivity  value  of  the  sea  surface  for  all  angles  [10]. 

Another  popular  but  realistic  method  of  computing  the  apparent  radiance  of  a  rough  sea  is  to 
create  a  surface  which  has  the  same  amplitudes  of  its  Fourier  components  as  a  rough  sea.  The 
method  is  called  the  power  spectral  approach  and  has  been  used  either  directly  or  indirectly  in  simu¬ 
lating  images  of  a  rough  sea  [11,12,13,3,14,15].  (See  [15]  for  a  description  of  the  technique.)  The 
power  spectral  method  is  computationally  intensive  since  a  large  number  of  facets  must  be  considered. 
Furthermore,  a  lot  of  surfaces  must  be  generated  to  compute  any  statistical  information  about  the  sea 
surface  in  question  for  each  wind  speed.  Analyzing  the  power  spectral  model  is  indeed  difficult.  Our 
probabilistic  approach,  however,  immediately  generates  the  statistics  of  the  rough  sea  surface  with  a 
minimum  amount  of  computational  effort.  The  physics  is  contained  in  the  model  at  a  glance,  thus 
making  it  easy  to  analyze  and  verify  in  certain  trivial  cases. 
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Finally,  our  slope  magnitude  model  of  the  sea  surface  is  a  two  dimensional  model  based  upon  on 
linear  superposition  of  waves.  Real  sea  surfaces,  however,  have  nonlinear  waves,  as  pointed  out  in 
116,171.  In  particular,  Huang  and  co-workers  [7,8]  have  derived  the  probability  slope  distribution  for 
a  nonlinear  wave  in  one  dimension.  In  [16],  a  small  parameter,  ok,  is  introduced  where  a  is  the  stan¬ 
dard  deviation  of  the  surface  elevation  and  k  is  the  wave  number.  Using  the  assumption  that  the 
parameter  is  small,  the  probability  slope  density  can  be  shown  to  be  normal  to  order  (ok)2.  There¬ 
fore,  the  generalization  to  a  nonlinear  sea  model  should  not  impact  greatly  on  our  results. 
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Appendix  A 

The  starting  point  is  to  consider  the  function  z=h(x,y)  defined  to  be  the  water  height  above 
mean  sea  level  at  the  position  (a:  ,y)  in  the  horizontal  plane.  The  function  h  may  be  regarded  to  be  a 
random  function;  that  is,  h(x,y)  is  to  be  regarded  as  a  random  variable  for  each  value  of  x  and  v. 
Corresponding  to  this  type  of  random  variable  is  a  probability  density  function.  Because  the  surface 
of  the  ocean  is  translationally  invariant  in  a  statistical  sense,  the  density  does  not  depend  on  x  and  v . 
Accordingly,  we  denote  it  by  ph.  The  value  of  h(x,y)  is  the  sum  of  contributions  from  waves  of  a 
large  range  of  wavelengths.  To  the  extent  that  the  contributions  of  the  different  wavelengths  are  sta¬ 
tistically  independent,  the  central  limit  theorem  applies  and  ph  is  of  normal  form.  This  assumption  of 
independent  contributions  is  an  approximation  which  breaks  down  for  large  values  of  h  due  to  non¬ 
linear  hydrodynamic  interactions.  However,  it  is  a  good  approximation  for  our  purposes  here.  The 
wave  slope  is  the  gradient  of  h .  The  statistics  of  spatial  derivatives  can  be  calculated  by  the  use  of 
finite  differences,  to  wit: 

h(x',y)  ~  h(x,y) 
dx  x'  -  y 

for  x'  sufficiently  close  to  x.  (A  similar  relation  holds  for  dhldy.)  It  has  been  established  that 
wavelengths  shorter  than  a  certain  cutoff  contribute  very  little  to  the  slope  statistics.  Consequently 
Eq.  (Al)  may  be  considered  valid  as  long  as  |jcy  — jc  (  is  less  than  the  cutoff  length.  If  the  anisotropy 
introduced  by  the  wind  direction  is  ignored,  h(x' ,y)  and  h(x,y)  are  identically  normally  distributed 
random  variables.  They  are  correlated  with  a  correlation  length  depending  on  |jc'-.r|.  Since  the 
difference  of  two  correlated  normally  distributed  random  variables  is  itself  a  normally  distributed  ran¬ 
dom  variable,  it  follows  that  dh/dx  and  dhldy  are  normally  distributed.  The  wave  slope,  s.  is 
defined  to  be  magnitude  of  the  gradient  of  h ;  i.e.. 


Since  dh/dx  and  dhldy  are  normal,  independent  with  zero  mean  and  equal  variance,  then  the  random 
variable  s  has  a  Rayleigh  density,  given  by 

ps(s)<xse~*'/r  (A3) 

where  <r  is  the  slope  variance.  It  is  more  convenient  to  work  with  a  description  in  terms  of  spherical 
coordinates  on  a  unit  sphere.  We  introduce  the  angle  coordinates.  6,  <t> ,  to  describe  the  direction  of 
the  unit  normal,  n ,  to  the  facet  face.  The  angle  9  is  the  angle  between  and  the  normal  to  the  sea  sur¬ 
face  (the  z  direction)  while  <j>  is  the  azimuthal  angle  defined  such  that  it  is  zero  when  n  is  in  the  c  _r 
plane.  Using  the  fact  that  s  -  tan#,  it  follows  from  a  change  of  variables  that 

p„(9,<t>)  =  (1/2t)2 cr  tan  9  seer  (A4) 

It  is  easy  to  see  how  the  nonisotropic  case  in  which  the  upwind  and  crosswind  variances  are 
unequal  generalizes  the  slope  probability  distribution.  Assume  that  the  .r  and  v  directions  correspond 
to  the  upwind  and  crosswind  directions,  respectively,  and  also  suppose  that  the  mean  slopes  are  zero 
in  both  x  and  y  directions.  A  normal  vector  to  the  surface  at  (x ,y )  is  given  by 

N(x.y)  =  (dh/dx.  dh/dyA)  =  [jr  +  l)l/2n(x ,y  >.  <A5) 
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Decomposing  the  vector  N(x,y)  into  spherical  coordinates  (0,<£)  on  the  sphere  of  radius  [5 2  +  1]1/2, 
it  follows  that  sin  0  =*  s /[s  2  +  1]1/2  and  tan  0  =  s . 
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Reasoning  as  before,  the  slope  density,  p(dh/dx,  dh  I  by),  is  normal  in  both  arguments  with 
standard  deviations  ox  ,ay ,  and  correlation  coefficient  r .  We  make  the  change  variables  by  letting 

dh  /dx  =  “^s2  +  1  sin  6  cos  <£  =  tan  6  cos  <t>  (A6) 

dh  !dy  =  ^ s 2  +  1  sin  6  sin  <t>  =  tan  $  sin  <t>. 

Since  the  Jacobian  of  the  transformation  is  sec2  6  tan  0,  the  new  density  in  terms  of  (0,<£)  is  given  by 

-  ...  , .  sec2  0  tan  0  j  tan2  0  "j 

Ai(0.d>)  =  - /  ,  exp  <  -  — - — E(<t>,r )  k  (A7) 

2Tajffyv(i  -  r2)  L  2d  "  r  )  J 

where 

r-/,  cos 2  <t>  ,  sin2  <t>  2r  cos  <*>  sin  <t>  _ 

c(<p,r) - —  +  - - -  (A8) 

(?x  &y  °x  °y 

and  n  is  Af/||AJ||.  Notice  that  if  ax  =  ay  and  r  =  0,  then  the  probability  density  is  given  in  terms 
of  0  only  and  reduces  to  Eq.  (A4). 


